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NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 
RESEARCH MEMORANDUM 


INVESTIGATION OF THE EFFECTS OF PROFILE SHAPE ON THE 
AERODYNAMIC AND STRUCTURAL CHARACTERISTICS OF = 
THIN, TWO-DIMENSIONAL AIRFOILS 
AT SUPERSONIC SPEEDS 


By Elliott D. Katzen, Donald M. Kuehn, 
and William A. Hill, Jr. 


SUMMARY 


In order to determine the effects of thickness, trailing-edge blunt- 
ness, boattailing, and forward profile on the aerodynamic characteristics 
of thin airfoils, and to provide a check on the available theoretical 
methods, 31 airfoils were tested. Тһе airfoils were 2, l, and 6 percent 
thick and were tested at Mach numbers of 1.45 and 1.98 at Reynolds numbers 
of 1.0, 2.0, and 3.5 million in a clean condition and at a Reynolds number 
of 3.5 million with transition fixed. 


The aerodynamic advantage of very thin airfoils was shown by a rapid 
increase of maximum lift-drag ratio (e.g., from 5.8 to 14.4 at M = 1.45 
and a Reynolds number of 1.0 million) as the airfoil thickness ratio was 
decreased from 6 to 2 percent. Increased trailing-edge bluntness of the 
6-percent-thick airfoils caused a small decrease in maximum lift-drag 
ratio but a large increase in section modulus; for the 2-percent-thick 
airfoils, increased bluntness caused a large decrease in maximum lift- 
drag ratio with only a small increase in section modulus. This has spe- 
cial significance for propeller designers in that it indicates that 
propellers whose blade elements operate at supersonic speeds should have 
blunt trailing edges for the thick sections near the hub and relatively 
sharp trailing edges for the thin sections near the tip. The importance 
of maintaining a laminar boundary layer on very thin airfoils was shown 
by the decrease of maximum lift-drag ratio from 14.4 to 11.9 caused by 
fixing transition at the leading edge of a sharp-trailing-edge 2-percent- 
thick airfoil at а Mach number of 1.45. The effects of the different 
forward profiles and changes in boatteailing were generally such that а 
reduction in profile area reduced the minimum drag coefficient. The 
center of pressure of the airfoils moved forward with increased thickness 
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ratio and moved aft with increased trailing-edge bluntness. Available 
theoretical methods were adequate for calculating the lift and pitching 
moment of the airfoils under all conditions of the tests. Тһе theoreti- 
eal methods for calculating foredrag and correlation curves for estimating 
base pressure were adequate for predicting the total drag when the tran- 
sition position was known. 


INTRODUCTION 


Although thin airfoils have inherent structural limitations, they 
are an aerodynamic necessity for wings, control surfaces, and propellers 
in order to attain practical lift-drag ratios for supersonic airplanes 
and missiles. The structural limitations of thin airfoils can be allevi- 
ated by making the trailing edges blunt and by increasing the area of the 
boattailed and forward parts of the profile. However, the aerodynamic 
penalties, if any, involved in these increases in profile area must be 
evaluated and balanced against the structural improvements they afford. 
The best means of assessing bluntness and other thickness effects, with- 
out having the results obscured by plan-form effects, is by the study of 
two-dimensional data. 


Most of the experimental data available on blunt-trailing-edge air- 
foils (e.g., refs. l through 5) pertain to thickness ratios greater than 
5 percent. Since gains in aerodynamic and propulsive efficiency are to 
be expected with further reductions in thickness ratio, data are required 
for smaller thickness ratios. For these very thin airfoils, the boundary- 
layer thickness is often the same order of magnitude as the airfoil thick- 
ness and large viscous effects are to be expected, particularly near the 
base of the model. Thus, experimental evaluation of the available invis- 
eid theories, such as the linear, second-order, and shock-expansion 
theories, and of the analytical results of references 1, 2, 6, 7, 8, and 
9 is required. | 


The purpose of the present investigation, therefore, is to provide 
experimental data on thin, blunt-trailing-edge airfoils for evaluating 
the effects of thickness, trailing-edge bluntness, boattailing, and for- 
ward profiles, &nd to furnish a check on the available theoretical 
results. To accomplish these aims, 31 airfoils were tested in the Ames 
l- by 3-foot supersonic wind tunnel No. 1 at Mach numbers of 1.45 and 
1.98 and at Reynolds numbers of 1.0 to 3.5 million. 


SYMBOLS : 7 
A cross-sectional area of airfoil profile, sq in. 
b Length of boattailed section, in. 
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airfoil chord, in. 
section drag coefficient 


section base-drag coefficient 
section minimum drag coefficient 
section minimum foredrag coefficient 


wing drag coefficient 
section Lift coefficient 


rate of change of lift coefficient with angle of attack at 
zero angle of attack, per radian 


section lift coefficient for maximum lift-drag ratio 


wing lift coefficient 

section pitching-moment coefficient about midchord position 
wing pitching-=-moment coefficient about midchord position 
trailing~edge thickness, in. 


section moment of inertia about chord line, in.* 


section modulus, in. š 


section maximum lift-drag ratio 


Mach number 


ratio of base pressure to free-stream static pressure 


Reynolds number based on chord 


maximum thickness of airfoil, in. 
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ΔΚ. distance of center of pressure forward of midchord position 
in chord lengths 

α angle of attack, deg 

B ΜΕ L1 

Y ratio of specific heats (1.400 for air) 


EXPERIMENTAL CONSIDERATIONS 


Apparatus 


The tests were conducted in the Ames l- by 3-foot supersonic wind 
tunnel No. 1 This single-return, eontinuous-operation, variable-pressure 
wind tunnel has a Mach number range of 1.2 to 2.5. The Mach number is 
changed by varying the wall contour by use of flexible plates which com- 
prise the top and bottom walls of the tunnel. 


The side-support balance used to measure the &erodynamic forces, and 
a typical two-dimensional model installation are shown in figure 1. Since 
this is the first report in which data obtained with this balance have 
appeared, a more complete description of the balance is given than would 
normally be the case. The two-dimensional installation includes a 
through-span model, two boundary-layer plates, and two complete balance 
units, one located on each side of the tunnel (only one balance unit is 
shown in figure 1). The model is supported at each end by the balance 
floating beam, which is, in turn, supported entirely by 8 set of six 
(the yawing-emoment gage is not shown) interchangeable ring-type strain 
gages. These gages are attached to the main body of the balance through 
insulated supports. The complete unit of test model, floating beam, and 
ring gages comes in contact with the main body of the balance only through 
these insuleted supports, thereby permitting electrical contact to detect 
fouling of the unit. The balance is rotated by an electric motor driving 
the worm gear mechanism on the main body of the balance. An airtight 
drum fits over the entire balance unit and is vented to free-stream static 
pressure. Calibration of the balance showed that interaction effects 
between components of the balance were negligible. 


Tunnel-wall boundary-layer effects are eliminated by the use of two 
parallel boundary-layer plates 8 inches apart which are alined with the 
tunnel free stream and which remain fixed with respect to the main tunnel 
walls. The spacing between the tunnel wall and the boundary-layer plate 
is such that all side-wall tunnel boundary layer flows behind the plate. 
Circular plates in the boundary-layer plates rotate with the model so as 
to maintain a uniform, preset clearance around the model surface and 
model shanks. The surfaces of the circular plate and the boundary-layer 
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plate remain flush to within 0.0015 inch. A fairing which completely 
envelopes the model shank, thereby relieving it of any air load, is an 
integral part of the rotating circular plate and is connected solidly to 
the balance housing. The model shank is pinned directly to the floating 
beam. Clearance around the model surface and model shank was maintained 
at all times so that the full air load on the wing would be transmitted 
to the ring gages. The clearance between the wing surface and the 
boundary-layer plate was preset at approximately 0.030 inch for all the 
two-dimensional wings (see fig. 2(а)). An auxiliary test arrangement, 
used only with а semispan wing, which eliminates the 0.030-inch gap by 
using a 0.007-inch-thick rubber seal is shown in figure 2(b). The semi- 
span wing uses only one balance and boundary-layer plate; the other 
boundary-layer plate is replaced by a conventional 18-inch-diameter 
window. 


Models 


All the airfoil sections used in this investigation are shown in 
figure 3. Each of the 31 sections was derived from one of four basic 
shapes (having three different forward profiles); namely, (1) biconvex 
(i.e., two circular arcs), t/e = 0.04, (2) biconvex to the midchord point 
with the maximum thickness constant to the trailing edge, t/c = 0.04, 

(3) biconvex to the one-thid chord point with the maximum thickness 
constant to the trailing edge, t/c = 0.02, 0.04, and 0.06, and (4) NACA 
16-00h. Seven individual sections were obtained from each of the four 
blunt airfoils by altering the basic profile, as shown in figure 3(b). 

In each case, the amounts of bluntness were 100, 60, 30, and O percent 

of the airfoil maximum thickness. The boattail parameter used in the 
investigation was the ratio of the length of boattail section to the 

chord length. For the basic airfoils biconvex to ο/3, the b/ ο ratios 

were 0.05 and 0.33; for the basic airfoil biconvex to c/2, the b/c ratios 
were 0.05 and 0.50. The biconvex airfoil was not altered. One alteration 
was made to the NACA 16-004 airfoil by increasing the trailing-edge blunt- 
ness to 30 percent with the boattail a straight line from the base tangent 
to the original profile. Each of these airfoils had a chord Length of 

6 inches and a span of 8 inches between boundary-layer plates. 


Two other wings were used to evaluate quantities necessary for the 
complete analysis of the main airfoil data. They were (1) a 3/ h-ineh- 
chord wedge with a base height of 0.180 inch and a span of 8 inches 
between boundary-layer plates, (2) a 6-inch-chord, aspect ratio 2, double- 
wedge, 6-percent-thick wing with a gap that could be sealed (see fig. 2). 


The airfoils with blunt trailing edges had a 0.020-inch-diameter 
base-pressure orifice located symmetrically in the base with the pressure 
lead brought out at the root of the wing. All leading edges, except for 
the NACA 16-004 airfoil, and all sharp trailing edges had radii of 
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approximately 0.002 inch for practical reasons of construction and main= 
tenance. The ordinates for the NACA 16-004 section are given in table I. 


Test Procedure 


The airfoils were tested at Mach numbers of 1.45 and 1.98 and at 
Reynolds numbers of 1.0, 2.0, and 3.5 million at each Mach number with 
the airfoil surface clean. At a Reynolds number of 3.5 million, the air- 
foils were also tested with a l/l-inch-wide salt band starting 1/4 inch 
from the leading edge. At a Mach number of 1.45, the angle-of-attack 
range was limited to +60 due to the movement of the diffuser normal shock 
wave into the test section of the wind tunnel. However, at & Mach number 
of 1.98 an angle of attack of +18° was attained before this phenomenon 
occurred. Simultaneous measurements of total force and base pressure 
were made at each angle of attack. For all tests the humidity of the 
tunnel air was held to а value of less than 0.0003 pound of water vapor 
per pound of dry air which is sufficient to reduce condensation effects 
to & negligible amount. 


Since a salt band was used to promote boundary-layer transition, 
it was necessary to evaluate the wave drag of the salt. For this purpose 
8, 3/l-inch-chord wedge (see fig. 3) was tested &t zero angle of attack 
at both Mach numbers with and without the salt band. An attempt was made 
to standardize the size of salt crystals used by sifting through two 
screens of different mesh size thus eliminating all the very large and the 
very small crystals. Several repeat runs were made to check the consist- 
ency with which the salt bands could be applied. The tests indicated 
that the salt-band results could be repeated within the values of uncer- 
tainty in drag coefficient stated later in the section "Accuracy of рафа"; 
that is, the differences between runs with the salt band were no larger 
than differences between repeat runs without the salt band. For tests 
of the wedge with and without salt, the total force and.the base pressure 
were measured and the results corrected to foredrag. The difference 
between the foredrag results with and without salt was taken as the cor- 
rection for the wave drag of the salt. | И 


The effect of the clearance around the airfoil at the wing root оп 
the aerodynamic forces was evaluated by using identical wing shapes with 
the gap around one of the wings sealed with а O.OO0T7-ineh-thick rubber 
seal. The change in strain-gage calibration caused by that portion of 
the applied force required to stretch the elastic rubber seal was deter- 
mined and found to be negligible by a static-force calibration with model 
and balance installed in the tunnel. The double-wedge, semispan wings 
used for this check were run at the same test conditions as were the two- 
dimensional airfoils. Results of the gap tests, described later, showed 
that the effect of the gap was insignificant. 
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Reduction of Data 


All force data were reduced to the usual coefficient form for drag, 
lift, and pitching moment. The reference length and area used for the 
coefficients were chord length and plan-form area, respectively. Та 
addition to the total drag for all wings, the foredrag has also been pre- 
Sented in some instances. Base-pressure coefficients are given in the 
form of the ratio of base pressure to the tunnel free-stream static 
pressure. 


Corrections to Experimental Results 


Corrections to the measured force coefficients were computed for the 
wave drag of the salt band, for irregularities of the free-stream pressure 
and stream angle, and for the change in angle of attack due to elastic 
deformation of model and balance. 


The correction (described in detail in the section "Test Procedure") 
which was applied to the main airfoil data for the wave drag of the salt 
band is an average of several runs made with the 3/4-inch-chord wedge. 
The resulting values, expressed in drag coefficient form, are 0.0035 at 
a Mach number of 1.45 and 0.0022 at a Mach number of 1.98. This correc- 
tion was determined at zero angle of attack only, but it was applied as 
а drag-force correction at all angles of attack. It is believed, however, 
that the error introduced by this method of correction is of no practical 
consequence. 


Since the airfoil was rotated about its midchord point and therefore 
always remained in nearly the same region of the tunnel stream, a single 
correction for stream angle sufficed throughout the angle-of-attack range. 
This stream-angle correction, which was small, and, in general, different 
for the lift and moment curves, included the amount that the model angle 
Of attack was in error at the initial installation and any asymmetries in 
construction. The correction was applied very simply by shifting the 
lift and moment curves so that they would pass through the origin of the 
axes. Corrections to the force and moment coefficients due to free-stream 
pressure gradients were negligible for the airfoils of such small thick- 
ness ratios as were used in this investigation. The combined elastic 
deformation of the model and balance for all wings was negligible for all 
test conditions. 
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Accuracy of Data 


The accuracy of the final parameters used in the data analysis has 
been estimated by considering the known accuracy of the individual quan- 
tities used in determining these final values. The total uncertainty is 
given by the square root of the sum of the squares of the individual 
uncertainties. These uncertainties are given as + increments in the 
following table: 


1.0 
3.2 
1.0 
3.2 
1.0 
3.3 
1.0 
3.5 


t.0020 





Repeatability of data was checked by making several runs with a given 
model and was found to be consistent with the above values. In addition 
to the above quantities, another source of uncertainty is introduced by 
the gap around the airfoil at the root of the wing. The effect of the 
gap on the aerodynamic characteristics was estimated from the semispan 
data shown in figure 4. Since the difference between the runs with and 
without gap was no more than that for a given model with a gap run twice, 
the effect of gap was negligible. 


THEORETICAL CONSIDERATIONS 
Lift and Pitching Moment 
Two theoretical methods of calculating the Lift and pitching moment 
were used in the present report. Second-order theory was used for cal- 
culating the slopes at а = О, and shock-expansion theory was used at the 
higher angles of attack where nonlinearities are to be expected and 


greater accuracy than that of second-order theory is required. The 
second-order formulas used for the slopes at а = О (ref. 1) were: 


de c 
3 2h 
de = 20 (a Ич 3) (2) 
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where 
Ci = = (2) 
Je - 1 
_ (7 + 1)M? - M(M2 - 1) 
Co zd 2(м2 Е 1)2 (3) 
апа 


2 
de πας (1) 


Formulas 1 and 4 neglect the base pressure but it was shown in reference 1 
that the error involved in excluding the base pressure was negligible for 
lift and pitching moment. 


Drag 


The drag of the airfoils was predicted by adding calculated values 
of the pressure foredrag, skin-friction drag, and the base drag. The 
pressure foredrag was calculated by shock-expansion theory for all the 
airfoils. For the two round-leading-edge sections (NACA 16-004 airfoils 
with h/t = О and h/t = 0.3), the calculations were made by approximating 
the leading edge with the largest wedge angle for which the leading-edge 
shock wave remained attached. Laminar skin-friction coefficients were 
estimated from the Blasius flat-plate incompressible theory since differ- 
ences are negligible at the test Mach numbers between this theory and the 
more accurate theories which account for compressibility. Turbulent skin- 
friction coefficients were estimated from Cope's theory (ref. 10) since 
various experimental results summarized by Chapman and Kester in reference 
11 have shown this theory to be as accurate as any of the available theo- 
ries at the Mach numbers of the test. No estimate was made of the skin- 
friction coefficients for the Reynolds numbers for which natural transition 
and part laminar, part turbulent boundary layers were indicated on the 
airfoils because the location of transition was not known. The base drag 
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Was estimated from the correlation plots of reference 12. Both the base 
drag and the skin-friction drag were assumed constant with angle of attack. 


Lift-Drag Ratio and Optimum Lift. Coefficient 


The maximum Lift-drag ratio was calculated from the following second- 
order equation of reference l which is applicable for small values of h/c: 


1/2 
ΕΊΣ; 
ON = a n Cic (5) 


The optimum lift coefficient was calculated from 


est = (203. πα. G + Ë as) (6) 


which was derived from the formulas of reference 1. 
RESULTS AND DISCUSSION 


Results of tests of the 31 airfoils are presented in table II in the 
form of lift, drag, and pitching-moment coefficients and base-pressure 
ratios &s & function of angle of attack. The coefficients sre averages 
of the values measured for both positive and negative angles since the 
airfoils were symmetrical about the chord plane and the differences for 
the positive and negative angles were small and consistent with the uncer- 
tainties in the data listed previously. Since the data are tabulated, 
the only basic plots shown (figs. 5 through &) are typical lift, pitching- 
moment, drag, and lift-drag-ratio curves. The trends of the data with 
the parameters of the test (t/c, h/t, b/e, forward profile, M and R) and 
theoretical and experimental correlation plots are shown in figures 9 
through 16. The experimental results of the test, with theoretical pre- 
dictions for comparison, are summarized in table III as lift-curve and . 
moment-curve slopes at a= 0, minimum drag coefficient, maximum lift-drag 


er 
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ratio, and lift coefficient for maximum lift-drag ratio. The results are 
discussed in the following paragraphs. 


Lift 


Typical lift curves are shown in figure 5. The figure shows that the 
lift at any given angle increased with trailing-edge thickness and was 
relatively unaffected by airfoil thickness ratio as predicted by theory. 
The predicted departure from linearity with angle of attack was slightly 
less than that measured. 


Examination of table III shows that the normalized lift-curve slopes, 
Вет, аге principally а function of h/c and that there was no consistent 


variation of Вс with t/e, b/c, forward profile, В and M. The nor- 
ba 


malized lift-curve slopes are plotted against h/c in figure 9 for all 
the models. The mean of the data closely followed the predicted increase 
in Вет with. h/c. The spread in the dat& due to variations in tfe 


and the other parameters mentioned above did not produce deviations from 
the predicted curve by more than +5 percent, with but few exceptions. 


Pitching Moment and Center of Pressure 


Typical pitching-moment coefficients are shown as a function of Lift 
coefficient in figure 6. The data are for sharp-trailing-edge airfoils 
2 and 6 percent thick in figure 6(a) and are for 6-percent-thick airfoils 
with both sharp and fully blunt trailing edges in figure 6(b). The figure 
shows that at a given lift coefficient and trailing-edge thickness ratio, 
the effect of increased airfoil thickness ratio was to increase the 
pitching-moment coefficient. However, for given thickness ratio, the 
pitching-moment coefficient decreased with increased trailing-edge thick- 
ness. Although the predicted pitching-moment curve was slightly more 
nonlinear with lift for the 6-percent than the 2-percent-thick airfoils 
and the reverse occurred experimentally, the agreement between theory and 
experiment was good. 


The foregoing trends of lift and pitching-moment coefficient are 
reflected in the forward movement of the center of pressure (fig. 10) 
with increased airfoil thickness ratio and in rearward movement with 
increased trailing-edge bluntness. For the sharp-trailing-edge airfoils 
shown, the center of pressure moved from 2 to 8 percent of the chord for- 
ward of the midchord position at M = 1.45 with an increase in airfoil 
thickness ratio from 2 to 6 percent. Illustrating the opposite trend 
with increased trailing-edge bluntness, the center of pressure shifted 
from 8 to 4 percent of the chord forward of the midchord position for the 
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6-percent-thick airfoil at M = 1.45 with an increase from zero to full 
bluntness. The center of pressure was generally 1 to 2 percent of the 
chord further forward at М = 1.45 than at М = 1.98. This was a trend 
not predicted by theory. The effect of b/e and of forward profile, 

with t/e and h/t fixed, was generally to shift the center of presgure 
forward with increased profile area. The variation of the center of pres- 
Sure with Reynolds number can be seen from table III to be less than +0.5 
percent of the chord in most cases. 


It can be seen from figure 11 that the predicted and experimental 
center-of-pressure positions usually agreed within 1 or 2 percent of the 
chord, with better agreement at М = 1.98 than at М = 1.45. 


An example of how two-dimensional section data can be combined with 
linear-theory tip effects to predict the pitching moment of a finite 
aspect ratio wing is shown in figure 4. It can be seen that superimposing 
thickness effects, as given by second-order theory, on linear-theory tip 
effects improves the prediction of the pitching moment over that which 
would be calculated with thickness effects neglected. 


Minimum Drag Coefficient 


The variation of minimum drag coefficient with airfoil thickness 
ratio and trailing-edge thickness ratio is shown in figure 12. Data 
included pertain to various forward-facing profiles and boattail condi- 
tions, with transition fixed at R=3.5xl0° and with the airfoils in a 
clean condition at R-1.0x109. The data manifest a rapid increase in 
Cain With increased t/c or h/t and a decrease in Canin With increased 


Mach number. The percentage increase in Canin With increased h/t was 
larger for the 2-percent than the 6-percent-thick airfoil, as would be 
inferred from the theoretical optimum-airfoil results of reference 5. At 
M = 1.45, Cü,,, increased approximately HO percent (or 0.012) for | 
6-percent-thick &irfoils and 100 percent (or 0.006) for 2-percent-thick 
airfoils, with an increase in h/t from О to 1.0. The effects of changes 
in the basic airfoil and in b/e were such that, generally speaking, 
reduction in profile area, regardless of position on the airfoil, reduced 


Cámin* 
The effect of fixing transition at the leading edge of the airfoils 
is also illustrated in figure 12. The increase in сач shown was 


usually due to an increase in base drag, as well as to an increase in 
skin friction. It can be seen from table III that the effect of increased 
Reynolds number on the airfoils in a clean condition was to increase 
Canin’ This indicates that the transition region was moving forward on 


the airfoil surface with increased Reynolds number because Ca | would 
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be expected to decrease with increased Reynolds number if the boundary 
layer had remained completely leminar. 


The theoretical and experimental Cain Were generally in good 


agreement, as can be seen in the correlation plots of figure 13. The air- 
foils for which the predicted drag coefficients were slightly higher than 
the experimental coefficients usually had sharp trailing edges (table III). 
This is probably a result of the shock wave near the trailing edge inter- 
acting with the boundary layer and causing increased pressures and reduced 
drag on the rear portion of the airfoils, as first pointed out in refer- 
ence 13. For the blunt-trailing-edge airfoils, the trailing-edge shock 
wave occurs further downstream, and the pressures on the rear portion of 
the airfoil apparently are not increased. The data of figure 13 are pre- 
sented as total-drag, foredrag, and b&se-drag coefficients for R-3.5x109 
with transition fixed and R-1.0x109 with the airfoils in a clean condi- 
tion. With transition fixed at R=3.5x10°, skin-friction and base-pressure 
coefficients corresponding to a turbulent boundary layer were used in com- 
puting the theoretical Canin’ With the airfoils in a clean condition at 


R=1.0x108, skin-friction and base-pressure coefficients corresponding to 
a laminar boundary layer were used in computing the theoretical Chain ° 


That the boundary layer was laminar over most of the airfoil is indicated 
by the good agreement between theory and experiment for the foredrag. Тһе 
increase in base drag over the predicted values suggests that transition 
is occurring near the base, thereby decreasing the experimental base- 
pressure ratio and increasing the base drag. Theoretical Camin was not 


calculated for R-2.0x108 and В=3.5х106 with the airfoils in a clean 
condition because the Location of transition on the airfoils was not known. 


Reference lO shows that base pressures at a given Mach number can 
be correlated on the basis of the ratio of the boundary-layer thickness 
at the base to the base height. Since the boundary-layer thickness at 
the base of the airfoils with transition fixed near the leading edge was 
undoubtedly different than for natural transition which did not occur 
near the leading edge, the base pressure with transition fixed would not 
be expected to coincide with the base pressure with natural transition. 
That they did not coincide can be seen from table II. At R-3.5x108 
the base-pressure ratios for natural transition (airfoils clean) were 
higher then those for transition fixed. However, the differences were 
not large and the differences were smaller at М = 1.98 than at M = 1.15. 
in general, figure 13 shows that errors in the base drag of thin airfoils 
do not seriously affect the prediction of the total drag. 
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Drag Due to Lift and Lift-Dreg Ratio — 


Typical variations of drag coefficient and lift-drag ratio with lift О E 
coefficient are shown in figures 7 and 8. The drag curves were parabolic = 
and the agreement between theory and experiment was good. The drag coef- 
ficient increased with increased trailing-edge bluntness at a given lift О Е 
coefficient; the increase occurred at а decreasing rate with increased Бо we 
lift coefficient. The lift-drag-ratio curves were more sharply peaked 2 
for the 2-percent than for the 6-percent-thick airfoils and for the sharp- ` | 
than for the fully blunt-trailing-edge airfoils. Furthermore, the lift Е 
coefficient for (2/a)max decreased with decreased t/c and h/t. These & 


trends naturally result from the behavior of Canin and are in good к ка 
agreement wlth theory. | de 


Maximum Lift-Drag Ratio = 


The variation of maximum lift-drag ratio with airfoil thickness xs 
ratio and trailing-edge thickness ratio is shown in figure lh. The aero- — . 
dynamic advantage of very thin airfoils is illustrated by the rapid - 
increase of ια as with decreased t/c for the sharp-trailing-edge 
airfoil in a clean condition at Б=1.0х10® and М = 1.45; the (17d) аз 
increased from 5.8 to 14.4 with a decrease in t/c from 6 to 2 percent. 
The increase of (1/4) „у With decreased t/c was larger for the sharp- 
then the blunt-trailing-edge airfoils, and the effect of increased h/t 
was to decrease λα). The decrease was small for the 6-percent-thick 
airfoils and became larger with decreased t/c. The effect of the dif- 
ferent forward profiles and changes in b/c were generally consistent 
with the Canin results in that reduction in profile area increased 


(1/4) ney? Increasing the Mach number from 1.45 to 1.98 caused a larger 
decrease in (l/d)max for the 2-percent than for the Ó6-percent-thick air- 
foils. It is also evident from figure lh that large gains in (2/4) ney 
can be achieved by maintaining a laminar boundary layer on very thin air- 


folls. For the 2-percent-thick sharp-trailing-edge airfoil at М = 1.45, 
(2/8) лах Wes 14.4 for the airfoils in a clean condition at R-1.0x108 


and 11.9 with transition fixed at R=3.5xl0%. ' > 


i il Í l.. 


In order that the various airfoils may be compared on the basis of " = 
both an serodynamic and a structural criterion, a plot of the variation 
of παλια. with section modulus is shown in figure 15. The figure 
shows that large amounts of bluntness of the 2-percent-thick airfoils 
resulted in large decreases in (1/ἄ)ῃο but only small increases in sec- 
tion modulus. Тп contrast, large increases in trailing-edge bluntness 
of the 6-percent-thick airfoils caused small decreases in (1/d)max and ë 
large inereases in section modulus. These trends would be essentially 


ee 
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the same if structural parameters other than section modulus, for instance, 
torsional rigidity or bending stiffness, had been chosen for comparison. 
Ihe results have special significance for propeller designers in that they 
indicate that propellers whose blade elements operate at supersonic speeds 
should have blunt trailing edges for the thick sections near the hub and 
relatively sharp trailing edges for the thin sections near the tip. It 
should be pointed out that airfoils designed to have the minimum wave drag 
for a given structural criterion (refs. 6 and 8) and for the Mach number 
and Reynolds number conditions of these tests have smaller leading-edge 
angles than the airfoils tested. For the 6-percent-thick airfoils espe- 
cially, it is believed that if airfoils with the same structural charac- 
teristics but with smaller leading-edge angles had been tested, higher 
maximum lift-drag ratios would have been obtained. 


The good agreement between theory and experiment for Cla and cq ; 
is mirrored in the good agreement for (12/4) тах shown in figure 16. 


CONCLUSIONS 


An investigation to provide experimental aerodynamic data at Mach 
numbers of 1.45 and 1.98 on thin, two-dimensional, blunt-trailing-edge 
airfoiis afforded the following conclusions: 


1. The aerodynamic advantage of very thin airfoils was shown by.a 
rapid increase of maximum lift-drag ratio with decreased airfoil thickness 
ratio. 


2. Increased trailing-edge bluntness of the Ó-percent-thick airfoils 
caused a small decrease in maximum lift-drag ratio and a large increase 
in section modulus; for the 2-percent-thick airfoils, increased bluntness 
caused a large decrease in maximum lift-drag ratio and а small increase 
in section modulus. 


3. The importance of maintaining a laminar boundary layer on very 
thin airfoils was shown by the decrease of maximum lift-drag ratio from 
14.4 to 11.9 caused by fixing transition at the leading edge of a sharp- 
trailing-edge 2-percent-thick airfoil at a Mach number of 1.45. 


4, The effects of different forward profiles and changes in boat- 
tailing were such that, generally, any reduction in profile area reduced 
the minimum drag coefficient. 


2. The center of pressure of the airfoils moved forward with 
increased thickness ratio and aft with increased trailing-edge bluntness. 


CUAL 3t 
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6. Available theoretical methods were adequate for calculating the 
lift and pitching moment of the airfoils under all conditions of the test. 
The theoretical methods for calculating foredrag and correlation curves 
for estimating base pressure were adequate for predicting the total drag 
when the position of transition was known. 


Ames Aeronautical Laboratory 
National Advisory Committee for Aeronautics 
Moffett Field, Calif., Feb. 8, 1954 
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TABLE I.- AIRFOIL ORDINATES 


Upper and lower surface 


МАСА 16-004 | NACA 16-00} 
Station . h/t=0 h/t-0.3 
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L.E. radius: 0.078 
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[Stations and ordinates given in percent of airfoil chord] 
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TABLE II.- LIFT, DRAG, AND PITCHING-MOMENT COEFFICIENTS AND 
BASE PRESSURE RATIO FOR AIRFOILS TESTED 
(a) Basic airfoil:  Biconvex to c/3, t/c = 0.02 
(1) h/t = 1.0, b/e = O 


R-3.5x10$ 5 R«2.0x109 Rel.0x10* 
Transition fixed Airfoil clean Airfoil clean 
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TABLE ΤΤ.- LIFT, DRAG, AND PITCHING-MOMENT COEFFICIENTS AND 
BASE PRESSURE RATIO FOR AIRFOILS TESTED - Continued 
(a) Basic airfoil: Biconvex to c/3, t/c = 0.02 - Continued 
(3) h/t = 0.3, b/e = 0.05 


R=3.5x105 R=3.5x10Š Re2.0x109 Rel.Ox109 
Transition rixed Airfoil clean Airfoil clean Airfoil clean 
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TABLE TI,- LIFT, DRAG, AND PITCHING-MOMENT COEFFICIENTS AND 
BASE PRESSURE RATIO FOR AIRFOILS TESTED - Continued 
(a) Basic airfoil:  Biconvex to с/3, t/e = 0.02 - Continued 
(5) h/t = 0.6, b/e = 0.33 


TUR 
a, Transition fixed Airfoll clean Airfoil clean Airfoil <a 
deg 
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TABLE II.- LIFT, DRAG, AND PITCHING-MOMENT COEFFICIENTS AND 
BASE PRESSURE RATIO FOR AIRFOILS TESTED ~ Continued 
(a) Basic airfoil:  Biconvex to c/3, t/c = 0.02 - Concluded 
(7) h/t = О, b/c = 0.33 


rm. | gum 
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(b) Basic airfoil: Biconvex to c/3, t/e = 0.04 
(1) h/t = 1.0, b/e = O 


R»3.5x108 - Rx3. 5x10* R-2.0x108 R=1.0x105 
Transition fixed Airfoil clean Airfoil clean : Airfoil clean 
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TABLE II.- LIFT, DRAG, AND PITCHING-MOMENT COEFFICIENTS AND 
BASE PRESSURE RATIO FOR AIRFOILS TESTED - Continued 
(b) Basic airfoil: Вісопуех to c/3, t/c = 0.04 - Continued 
(2) h/t = 0.6, b/e = 0.05 


R23. 5x105 R«3.5x108 R=2.0x10 R=1.0x10 
Transition fixed Airfoil clean Airfoil clean Airfoil clean 
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TABLE II.- LIFT, DRAG, AND PITCHING-MOMENT COEFFICIENTS AND 
BASE PRESSURE RATIO FOR AIRFOILS TESTED - Continued 
(b) Basic airfoil: Biconvex to c/3, t/c = 0.04 - Continued 
(4) n/t = 0, b/e = 0.05 


R»3.5x108 R=3.5x1o К=р.0х106 А-1 .0х10@ 
Transition fixed Airfoil. clean Airfoil clean Airfoil clean . 
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TABLE ΤΤ.- LIFT, DRAG, AND PITCHING-MOMENT COEFFICIENTS AND 
BASE PRESSURE RATIO FOR AIRFOILS TESTED - Continued 
(b) Basic airfoil: Biconvex to c/3, t/e = 0.04 - Concluded 
(6) h/t = 0.3, b/e = 0.33 


R=3.5x105 R=3.5x108 R=2.0x10° R=1.0x10 
Transition fixed Airfoil clean Airfoil clean Airfoil clean 
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TABLE Il.- LIFT, DRAG, AND PITCHING-MOMENT COEFFICIENTS AND 
BASE PRESSURE RATIO FOR AIRFOILS TESTED - Continued 
(c) Basic airfoil: Biconvex to c/3, t/c = 0.06 
(1) h/t = 1.0, b/e = 0 


Re}. 5x108 Re3, 5x10? R=2.0x10° В-1.0х108 
Transition fixed Airfoil clean Airfoil clean Airfoil clean 
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TABLE II.- LIFT, DRAG, AND PITCHING-MOMENT COEPFICIENTS AND 
BASE PRESSURE RATIO FOR AIRFOILS TESTED - Continued 
(c) Basic airfoil: Вісопуех te c/3, t/e = 0.06 - Continued 
(3) h/t = 0.3, b/e = 0.05 


Β.3. 5x10 Re3. 5x10 R-2.0x108 Re1.0x109 
Transition fixed Airfoil clean Airfoil clean Airfoil clean 
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TABLE II.- LIFT, DRAG, AND PITCHING-MOMENT COEFFICIENTS AND 
BASE PRESSURE RATIO FOR ATRFOILS TESTED - Continued 
(c) Basic airfoil: Biconvex to с/3, t/c = 0.06 ~ Continued 
(5) h/t = 0.6, b/c = 0.33 


Re3. 5x10? R=3. 5x10 В-2.охлоё R=1.0x10° 
Transition fixed Airfoil clean Airfoil clean Airfoil clean 
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TABLE II.- LIFT, DRAG, AND PITCHING-MOMENT COEFFICIENTS AND 
BASE PRESSURE RATIO FOR AIRFOILS TESTED - Continued 
(c) Basic airfoil: Biconvex to c/3, t/c = 0.06 - Concluded 
(7) h/t = O, Ъ/е = 0.33 


R»3. 5x108 R=3.5x10° R-2.0x108 
Transition fixed Airfoll clean Airfoil clean 
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airfoil: Biconvex to c/2, 
(1) h/t = 1.0, b/e = 0 
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TABLE II.- LIFT, DRAG, AND PITCHING-MOMENT COEFFICIENTS AND 
BASE PRESSURE RATIO FOR AIRFOILS TESTED - Continued 
(4) Basic airfoil: Biconvex to c/2, t/e = 0.04 - Continued 
(2) h/t = 0.6, b/e = 0.05 


R=3.5x10° R=3. 5x109 R-2.0x10° R=1.0x10° 
Transition fixed Airfoil clean- Airfoil clean Airfoil clean 
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TABLE II.- LIFT, DRAG, AND PITCHING-MOMENT CORFFICIENTS AND 

BASE PRESSURE RATIO FOR AIRFOILS TESTED - Continued | 

(d) Basic airfoil: Вісоруех to c/2, t/e = 0.04 - Continued 
(4) h/t = 0, b/e = 0.05 
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TABLE IT.- LIFT, DRAG, AND PITCHING-MOMENT COEFFICIENTS AND 
BASE PRESSURE RATIO FOR AIRFOILS TESTED - Continued 
(d) Basic airfoil: Biconvex to c/2, t/c = 0.04 - Concluded 
(6) h/t = ο. 3 b/e - 0.50 
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Trangition fixed Airfoil clean Airfoil clean Airfoil clean 
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TABLE II.- LIFT, DRAG, AND PITCHING-MOMENT CORFFICIENTS AND 
BASE PRESSURE RATIO FOR AIRFOILS TESTED - Continued 
(е) Basic airfoil: Вісопуех, t/c = 0.04 


R23.5x108 R23.5x109 R-2.0x109 R-1.0x109 
Transition fixed Airfoil clean Airfoil clean Airfoil clean 
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TABLE II.- LIFT, DRAG, AND PITCHING-MOMENT COEFFICIENTS AND 
BASE PRESSURE RATIO FOR ATRFOILS TESTED - Concluded 
(f) Basic airfoil: МАСА 16-004 - Concluded 
(2) h/t = 0.3 


R=3. 5x108 R=3.5X10° R-2.0x109 Rel.Ox108 
Transition fixed Airfoil clean Airfoil clean Airfoil clean 
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TABLE III.- SUMMARY OF RESULTS 






Basic airfotl: 
Configur- biconvex to c/3 
ou er 0.02 


h/t = 1.0, b/c = 0 
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TABLE III.- SUMMARY OF RESULTS = Continued 


Basic airfoil: 
biconvex to c/3 
t/o = 0.04 


h/t = 1.0, b/e = 0 
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TABLE III.- SUMMARY OF RESULTS = Continued 
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TABLE III.- SUMMARY OF RESULTS - Continued 


Basic airfoil: 
biconvex to с/о 
t/c = 0.0k 
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TABLE ΤΤΤ.- SUMMARY OF RESULTS - Concluded 


Вісоптех, t/e = 6.08, 
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Le~ Floating beam 
2.- Normal-force gage 


3," Rolling-moment gage 9.- Worm-gear drive mechanism 

h.- Pitching-moment gage 10.- Balance housing 

5e- Side-force gage ll.- Two-dimensional model 

6.- Pin connecting model 12.- Rotating circular plate 
Shank to floating beam 13.- Fairing 

Te- Сһога-Рогсе gage 14.- Boundary-layer plates 

8,- Model shank 15.- Main tunnel walls 


Figure l.- Side-support balance with a typical two-dimensional model 
installation. 
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(b) Model with seal, 





Figure 2.- Details of the junction between the wing and the boundary- 
layer plate. | | 
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b= length of boattail section 
h = trailing-edge thickness 
! = airfoil thickness 


(1) Biconvex 
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(2) Biconvex to c/2 
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(4) NACA 16-004 
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(6) Semispon model. 
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t/c: 06, AR = 2, c= 6” 


Figure 3.- Method of designating airfoils amd summary sketch of models 
(upper half shown, not to scale). 
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(a) Effect of gap on lift and pitching moment. 
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Figure 4,- Effect of gap on lift, drag, and pitching moment of semispan model; М = 1.98. 
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(b) Effect of gap on drag. 


Figure 4.- Concluded. 
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Basic airfoil: 
Biconvex to c/3 
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ema Shock-expansion theory 
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Typical variations of lift coefficient with angle of attack; 
М = 1.98, transition fixed, R-3.5x109. 
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(a) Thickness effect; h/t = 0, b/e = 0.33.(b) Treiling-edge bluntness effect; t/e = 0.06. 


M = 1.98, 


coefficient with lift coefficient; 
3.5х108. 


R= 


transition fixed, 


Figure 6.- Typical variations of pitching-moment 
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Basic airfoil: Biconvex to c/3 || fal | 
о 6.02, Mt х1, 5/6 = 0 
© Lt-.02, ht, b/=.33 үр δ]. 
E xx .06, h/t =I, 5/6 = О 
Т = 06 ҺО ὁ-.58 Д ЖЕ 
Filled symbols denote 
shock-expansion theory 








Drag coefficient, c, 
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Figure 7.- Typical drag curves; M = 1.98, transition fixed, R-3.5x108. 


ο ο | Eq Basic airfoil: Biconvex to c/3 
© 102, hZf=/ b/ =O 

MEN ΜΡ © 0702, h/! «0 be ».33 
iz s B Jc-06 ht, б/с” O 

Е + c=O6, h/t =O, ме = 33 


° E . 
||| | 1 ΠΡΙ. | геа ото denote 


















3 Ú ttt || pats | ESA shock- expansion theory 

Fae | 

š 

: ИШИ ЫС Ж 
ar μι LT ТЕ a | _ 





= 
wA ΓΓΗΓΓῚ uz LLL 
ТЕА i 


: [^ 
Lift coefficient, с; 


Figure 8.- Typical lift-drag ratio curves; M = 1.98, transition fixed, 
Rz3.5x108. 
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Figure 9.- Variation of lift-curve slope with trailing-edge bluntness for 
all models tested; R = 1.0, 2.0, and 3.5 million, M = 1.45 and 1.98. 
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Ф AA=LO A Biconvex to с/2, М6 ».04, b/c».50 
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(b) M = 1.98 


Figure 10.- Variation of center-of-pressure location with thickness ratio 
and tralling-edge thickness; transition fixed, R-3.5x108. 
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Figure 11.- Correlation between theoretical and experimental center-of- 
pressure locations for all the airfoils and Reynolds numbers. 
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Basie airfoil: Biconvex to с/3 —oO— Biconvex to c/3, Ме = 02, b/£- 33 
—Qo— 450, b/ç -.33 σ ме =.02, b/c=.05 
—0— 20:0, b/c =.05 ~—f}-— WO = OF b/t-.33 
------ 4416, b/e ».33 —o— We=.06, b/c-.33 
—o— 0/0510 b/c* Q σ МС =.06, δ/.-.05 
Filled symbols denote —4&— Biconvex to с/2, М6 =.04, Б =.50 
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Figure 12.- Variation of minimum drag coefficient with airfoil thickness 
ratio and trailing-edge thickness ratio; transition fixed, R=3.5x10® 
unless noted otherwise. 
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(а) M = 1.55, transition fixed, (b) М = 1.98, transition fixed, 
R-3.5x109. R-3.5x105. 


Figure 13.- Correlation between theoretical and experimental minimum 
drag coefficients. 
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(c) M = 1.45, model clean, (d) M = 1.98, model clean, 
R-1.0x109 ° R=1.0x10°. 


Figure 13.- Concluded. 
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Basic airfoil: Biconvex fo C/3 ____о Biconvex to c/3, М *.02, bfe ».33 
-----6 4-0 2/6 =.3.3 σ М =.02, МС =.05 
б h/-0 bc =.05 — — 8 Ме =.04 b/csz.33 
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---- h/!-LO, e MC*.06, b/cs.05 
Filled symbols indicate —.— А Biconvex to c/2, Мет 04 b/c».5 
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Figure 1l.- Variation of maximum lift-drag ratio and airfoil section 
modulus with thickness ratio and trailing-edge thickness; trsnsi- 
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Figure 15.- Variation of airfoil maximum lift-drag ratio with section 
modulus, 
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(b) M = 1.98, transition fixed, R-3.5x108. 


Figure 15.- Continued. 
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Figure 15.- Continued. 


PER СЕ 


Οἱ 


58 


Maximum lift-drag ratio, (Id) ox 


ADDL 


Pes 0g e— — М6 = OF cx /©=.06 





/3 


/2 





// 


10 


ο 02 OF 26 OF 


JO 
Section modulus, IAt/2) 


NACA RM A54B08a 


Basic airfoil 


© 
[3 
Φ 
Δ 
v 
a 


Biconvex fo с/3, МЕ =.02 
t/¢*.04 
c =.06 
Biconvex fo c/2, МС «.04 
Biconvex, c» OF 
МАСА 16-004 


Je {4 “6 


(a) М = 1.98, airfoils clean, В=1.0х108. : 


Figure 15.- Concluded. 
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Figure 16.- Correlation between theoretical and experimental maximum 
Lifte-drag ratios. 
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